The Olney transit corridor runs along Olney Ave between Chelten Avenue and Front Street. The 18 and 26 are the two routes that run along the corridor, both crossing North Broad Street past the Olney Transportation Center. The Transportation Center has a stop for the Broad Street Line, while there is also a stop for the Fox Station Regional Rail line at the end of the corridor adjacent to Olney & Front Street. There are length(subcorridors[[3]][[1]]) total bus stops along the corridor, with the average bus trip running past round(mean(subcorridors_dat[[4]][[1]]$n_stops)) stops spanning round(mean(subcorridors_dat[[4]][[1]]$distance_traveled),digits=2) miles.
# map of routes that exist in this corridor
leaflet() %>%
setView(lng = -75.14511, lat = 40.03905, zoom = 12) %>%
addProviderTiles(providers$Stamen.Toner) #%>%
#addPolylines(data = routes_w_ridership, color = "#4377bc", weight = 4, layerId = link_stop_data$fromto, opacity = 0.5)
# also add chart/table of global averages for context of subcorridors
These charts illustrate characteristics for the Olney corridor as a whole, such as average speed and ridership. Average speed, as well as riderhship, vary both by route and by time of day.
table_1 <- route_analytics %>% bind_rows(analytics %>% mutate(route_id = "Total"))
kable(table_1, booktabs = TRUE, align = 'c',format.args = list(big.mark = ","),digits=1) %>%
kable_styling(latex_options = "scale_down") %>%
row_spec(dim(table_1)[1], bold = T) %>% # format last row
column_spec(1, italic = T) %>% # format first column
scroll_box(width = "100%", height = "300px")
| route_id | daily_ridership | trips | routes_served | service_hours | riders_per_hour | on_off | dwell_observed_mean | dwell_predicted_mean | dwell_hybrid_mean | dwell_per_onoff | onoff_per_trip | onoff_per_tripstop | avg_segment_speed | avg_speed_10_pct | avg_speed_25_pct | avg_speed_75_pct | avg_speed_90_pct |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 18 | 10,592 | 235 | 18 | 30.1 | 351.6 | 11,448.4 | 0 | 216.0 | 216.0 | 0.2 | 48.7 | 2.1 | 11.2 | 9.1 | 9.6 | 11.9 | 14.2 |
| 26 | 8,414 | 229 | 26 | 55.5 | 151.5 | 9,839.4 | 0 | 199.3 | 199.3 | 0.2 | 43.0 | 1.8 | 10.9 | 9.4 | 9.9 | 11.6 | 13.2 |
| Total | 19,005 | 464 | 18, 26 | 85.7 | 221.9 | 21,287.8 | 0 | 207.8 | 207.8 | 0.2 | 45.9 | 1.9 | 11.0 | 9.2 | 9.8 | 11.7 | 14.0 |
# ggplot(route_analytics, aes(x = route_id, y = avg_segment_speed)) +
# geom_bar(stat = "identity", fill="skyblue", alpha=0.7) +
# geom_errorbar(data = route_analytics, stat = "identity", ymin = route_analytics$avg_speed_10_pct, ymax = route_analytics$avg_speed_90_pct, colour="orange", alpha=0.6, size=1.3, width = 0.4) +
# scale_y_continuous(name = "Speed (MPH) - Includes Dwell Time", n.breaks = 8, limits = c(0, max(route_analytics$avg_speed_90_pct) * 1.1)) +
# scale_x_discrete(name = "Route Number") +
# labs(title = paste0("The Average Bus Travels at ", round(mean(analytics$avg_segment_speed), 1), " MPH on the Corridor")) +
# theme(text = element_text(size = 9))
plot_hourly_speed(subcorridors_dat,"OL")
ggplot(binned_analytics, aes(x = timeframe, y = avg_segment_speed)) +
geom_bar(stat = "identity", fill="skyblue", alpha=0.7) +
geom_errorbar(data = binned_analytics, stat = "identity", ymin = binned_analytics$avg_speed_10_pct, ymax = binned_analytics$avg_speed_90_pct, colour="orange", alpha=0.6, size=1.3, width = 0.4) +
scale_y_continuous(name = "Speed (MPH) - Includes Dwell Time", n.breaks = 8, limits = c(0, max(binned_analytics$avg_speed_90_pct) * 1.1)) +
scale_x_discrete(name = element_blank()) +
labs(title = "Speed by Period (All Routes)")+
theme(text = element_text(size = 9))
ggplot(binned_analytics, aes(x = timeframe, y = daily_ridership, label = daily_ridership)) +
geom_bar(stat = "identity", fill="skyblue", alpha=0.7) +
geom_text(position=position_dodge(width=0.9), vjust = -0.75, size = 3) +
scale_y_continuous(name = "Total Riders Served", n.breaks = 5) + #, limits = c(0, max(binned_analytics$avg_speed_90_pct) * 1.1)) +
scale_x_discrete(name = element_blank()) +
labs(title = "Ridership by Period (All Routes)")+
theme(text = element_text(size = 9))
Here, we separate overall corridor statistics by route to better understand the makeup of Olney bus traffic. The first three graphs, which cover ridership (both average daily ridership and average ridership per hour), the number of trips, and number of service hours per route, aggregate for both directions of each route. Next, average hourly ridership is again shown, but divided by whether the bus was Eastbound or Westbound (as Olney runs E-W). Average corridor running time, or the time it takes for a bus to make one trip along the corridor, and average hourly speed are also differentiated by direction.
grid.arrange(ncol=2,
ggplot(route_analytics, aes(y = daily_ridership, x = route_id, label = daily_ridership)) +
geom_bar(stat = "identity", fill = "skyblue") +
geom_text(position=position_dodge(width=0.9), vjust = -0.75, size = 3) +
scale_y_continuous(name = "Total Riders Served", n.breaks = 8, limits = c(0, max(route_analytics$daily_ridership) * 1.1)) +
scale_x_discrete(name = "Route Number") +
labs(title = paste0("The Corridor Serves ", round(sum(analytics$daily_ridership)), " Riders per Day")) +
theme(text = element_text(size = 9)),
ggplot(route_analytics, aes(x=route_id, y = trips, label = trips)) +
geom_bar(position="dodge", stat="identity") +
geom_text(position=position_dodge(width=0.9), vjust = -0.75, size = 3) +
ylab("Daily Trips (2019)") +
scale_x_discrete(name = "Route Number") +
scale_fill_phl(palette = "main", discrete = T) +
labs(fill= "", title = paste0("The Corridor Serves ", sum(route_analytics$trips, na.rm = TRUE), " Trips per Day")) +
#theme_phl() +
theme(legend.position = "top")+
theme(text = element_text(size = 9))
)
ggplot(route_analytics, aes(x=route_id, y = service_hours, label = round(service_hours))) +
geom_bar(position="dodge", stat="identity") +
geom_text(position = "dodge", stat = "identity", vjust = -0.25) +
#ggtitle("Market and JFK serve over 1,100 Bus Trips per Day") +
ylab("Service Hours Per Day") +
scale_x_discrete(name = "Route Number") +
scale_fill_phl(palette = "main", discrete = T) +
labs(fill= "", title = paste0("The Corridor Serves ", round(sum(route_analytics$service_hours, na.rm = TRUE), 1), " Service Hours per Day")) +
#theme_phl() +
theme(legend.position = "top")+
theme(text = element_text(size = 9))
ggplot(hourly_route_analytics, aes(x = trip_hour, y = daily_ridership, group=route_id, fill=route_id)) +
geom_bar(position = "stack", stat = "identity") +
#geom_point(shape=21, size=3) +
scale_color_viridis(discrete = TRUE) +
scale_y_continuous(name = "Riders Served", n.breaks = 8) +
scale_x_continuous(name = "Hour", breaks = c(0:23)) +
labs(title = "Hourly Ridership by Route (all directions)") +
theme(text = element_text(size = 9))
plot_ridership_by_route_dir(subcorridors_dat, "OL")
ggplot(hourly_route_direction_analytics %>% subset(avg_run >=0), aes(x = trip_hour, y = avg_run, fill=direction_id)) +
geom_bar(position = "stack", stat = "identity") +
#geom_point(size=3) +
#scale_fill_viridis(discrete = TRUE) +
ylab("Average Running Time (Includes Dwell)") +
scale_x_continuous(name = "Hour", breaks = c(0:23)) +
labs(title = "Average End-to-End Running Time by Route/Direction") +
theme(text = element_text(size = 9)) +
facet_wrap(~route_id + direction_id, ncol = 2) +
#theme_phl(base_size = 9)
theme(text = element_text(size = 9), axis.text.x = element_text(angle = 90))
plot_speed_by_route_dir(subcorridors_dat, "OL")
Note: average speed and running time are calculated for the entire running time, including dwell times.
There are three main “sub-corridors” within the Olney corridor: Chelten to Broad, which is the section west of Broad St, Broad to 7th, and 7th to Front, which are both east of Broad St.
#stop_cb <- c("SEPTA372","SEPTA15911","SEPTA15794", "SEPTA15912", "SEPTA15789", "SEPTA15915", "SEPTA15586", "SEPTA15793", "SEPTA15913","SEPTA15791", "SEPTA15914", "SEPTA15792", "SEPTA15786", "SEPTA15916", "SEPTA15782", "SEPTA15917", "SEPTA15587", "SEPTA16979", "SEPTA15779", "SEPTA15918", "SEPTA15919")
plot_daily_ridership(subcorridors_dat, "CB")
plot_hourly_speed(subcorridors_dat, "CB")
plot_speed_by_route_dir(subcorridors_dat,"CB")
stop_b7 <- c("SEPTA15796", "SEPTA15908", "SEPTA15798", "SEPTA15907", "SEPTA15799", "SEPTA15800", "SEPTA16966", "SEPTA15814", "SEPTA16965", "SEPTA381", "SEPTA15795", "SEPTA373", "SEPTA15910")
trip_dat_b7 <- (find_trip_dat_v2(apc_data, stop_b7))
## [1] "Running passenger data for 464 trips"
analytics_b7 <- analyze_segment(trip_dat_b7)
stop_7f <- c("SEPTA16963", "SEPTA15815", "SEPTA16964", "SEPTA15817","SEPTA15819", "SEPTA16961", "SEPTA32312", "SEPTA16960", "SEPTA15820", "SEPTA16959", "SEPTA15822", "SEPTA15697", "SEPTA15899", "SEPTA15818", "SEPTA16962")
trip_dat_7f <- (find_trip_dat_v2(apc_data, stop_7f))
## [1] "Running passenger data for 452 trips"
analytics_7f <- analyze_segment(trip_dat_7f)